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The environment of high-redshift galaxies is characterized by both wind-driven out- 
flowing gas and gravitationally infalling streams. To investigate such galaxy-IGM in- 
terplay we have generated synthetic optical absorption line spectra piercing the volume 
surrounding a starbursting analog of a Lyman Break Galaxy selected in a z w 3 output 
from a SPH simulation, including a detailed treatment of mechanical feedback from 
winds. Distributions for several observable species (Hi, Cm, Civ, Sin, Sim, Siiv, Ovi, 
Ovn, and Ovm) have been derived by post-processing the simulation outputs. The hot 
wind material is characterized by the presence of high-ionization species such as Ovi, 
Ovn, and Ovm (the latter two observable only in X-ray bands); the colder (T < 10 5 5 
K) infalling streams can be instead identified by the combined presence of Sin, Sim, 
and Cm optical absorption together with Ovi that surrounds the cooler gas clumps. 
However, both line profile and Pixel Optical Depth analysis of the synthetic spectra 
show that the intergalactic filament in which the wind-blowing galaxy is embedded 
produces absorption signatures that closely mimic those of the wind environment. We 
conclude that it may be difficult to clearly identify wind-blowing galaxies and their 
complex gaseous environment at high redshift in optical QSO absorption-line spectra 
based solely on the observed ion absorption patterns. 
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1 INTRODUCTION 

Shortly after the first detection of quasi-stellar objects 
(quasars or QSOs), the presence of many narrow absorp- 
tion lines in the spe ctra of t h ese ob jects was recognized for 
the first time (e.g., iBahcalll (11966m ). It soon became clear 
that these lines ( "Lya forest" ) are produced by absorption 
of neutral hydrogen in gaseous structures that fill the inter- 
galactic space (the InterGalactic Medium, IGM). Today it is 
widely accepted that the neutral hydrogen Lya forest traces 
baryon density fluctuations associated with gravitational in- 
stability as part of the large-scale structure formation in the 
Universe. 

With the advent of more advanced spectrographs and 
larger telescopes, higher quality spectra were obtained 
and with great surprise it was discovered that the hy- 
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drogen forest is complemented by absorption lines from 
heavy elemen t s (ICowie et al. 1995). Sub sequent analysi s 
IIEllison et al. | i| 19991 ) . ISchave et al.l (|200ol ). ISongailal i|200ll) . 
ISimcod ( 20061 )) have shown that the IGM is enriched with 
metals at virtually any overdensity and redshift probed, as 
shown by the observation of resonance transitions associated 
to ions like Cm, Civ, Sin, Sim, Siiv, Ovi, and others; the 
enrichment of und erdense regions is still matter of debate 
|Aracil et al l |2004| L A major problem in our current un- 
derstanding of the IGM is to determine when and by what 
transport means metals have been able to travel to large 
distances from their productions sites (presumably, stars in 
galaxies. It is then clear that the answers to these questions 
have to come from the study of the physical galaxy-IGM 
interplay. 

One possible mechanism is dynamical stripping. During 
a close encounter of two or more galaxies, tidal forces can 
strip away from them part of the enriched gas, which is then 
dispersed into the IGM. This scenario is particularly impor- 
tant at high redshift, when the encounter rate is enhanced. 
However, this mechanism turns out to be too inefficient in 
reproducing the observed metal distribution with Hi col- 



2 A. P. M. Fangano et al. 



umn densities (|Aguirre et al.l l2001). A second mechanism is 
based on the idea that the energy /momentum deposited by 
multiple superno vae can drive an out fl ow of material from 
the h os t galaxy dDekel fc Silk! dl986h: iMac Low fc Ferraral 
1 19991) . IFerrara fc Tolstovl d2000h: iMadau, Ferrara fc ReeT 
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Dave fc Qppenheimerl (j2006l )). In this scenario, multi- 



ple supernova-driven bubbles merge together to form a 
"superbubble" ; the entrapped material would blow-out 
from the galaxy and vent the hot interior (largely enriched 
by freshly produced heavy elements) into the IGM giving 
rise to a large-scale outflow. Yet, the epoch when most 
of enric hment too k plac e is still unclea r. Stud ie s carried 
out b y ISongailal il200ll), ISchave et al] (120031 ), ISongailal 
|2005l ). iRvan- Weber. Pettini fc Madaul l |2006j ) have shown 
that the observed Civ absorption characteristics do not 
evolve significantly in the redshift range 2 < z < 6. This is 
surprising in view of the fact that most of the star-f ormation 
activit y in the universe took place at z < 5 (see iHopkinsl 
(|2004 ) for an updated compilation). The findings were 
interpreted as the results of an early [z > 6) metal pollution 
(often referred to as pre- enrichmen t) by the first generation 
(Pop HI) of stars and galaxie s l|Madau. Ferrara fc Reesl 
|200ll ). IScannapieco et aU l |2002j )). 

With the combined analysis of QSO and galaxy spec- 
tra it is possible to connect the properties of the in- 
tergalactic absorption lines with that of the for eg round 
gala xies. Using this techn ique lAdelberger et al.1 (|2003l ) 
and lAdelberger et ail (|2005l ) demonstrated that almost all 
galaxies at z > 2 can be associated with metal ab- 
sor bers, in agreement with the e arlier theoretical study 
by lHaehnelt. Steinmetz fc Rauchl l|l996[). and sub s equen t 
analysis by Desiacques et al. (l2004ir Mase lli et all " " 1 : 
iBertone fc White! (|2005l ^l . and iKawata fc Rauch 



2004), 



2007 



While this points directly to the hypothesis that those galax- 
ies are experiencing large scale outflows and hence polluting 
the IGM in situ, the actual mass loss rate and, most impor- 
tantly, whether they are able to escape the gravitational pul l 
of the galaxy are very hard to assess. ISimcoe et ail (|2006l ) 
have recently analyzed a line of-sight toward a quasar at z ~ 
2.7 and searched for a correlation between absorbers and 
nearby galaxies. By constructing detailed ionization mod- 
els for the different absorbers, they find that their observa- 
tions indicate a metal-enriched volume of ~ 100 — 200/1^ 
physical kpc around the analyzed galaxies. Furthermore, the 
sheet-like geometry of the absorbing clouds (AL < 1 kpc) 
as well as their high metallicity (Z > 0.1 — 0.3Zq) suggest 
an outflow scenario for their origi n . Con centrating on low 
redshifts, (z ~ 0.5) iBouche et alj l|2006l ). starting from a 
sample of 1806 Mgn absorbers, show that their equivalent 
widths (W > 2 A) are anticorrelated with the absorber halo- 
mass. Since this indicates that absorbing clouds are not in 
virial equilibrium with the host-galaxy, they concluded that 
these absorbers are most probably generated by a starburst 
activity and expelled via supernova-driven winds. Note that 
an early-epoch enrichment model is not excluded by this 
scenario since metals ejected by the first galaxies (z ~ 6-12) 
should also end up in ove r dense regions at lo wer redshifts 
l|Porciani fc Madaul (120051 ) . IScannapieco! (|2005h ). 

The complexity of absorbing cloud populations asks for 
a careful analysis and interpretation of observational data. 



For this aim, simulations can be quite helpful, especially 
in trying to identify the various gas components that con- 
tribute to the gas flow in the intergalactic environment of 
galaxies. However, simulations are not yet able to follow 
in detail the full range of physical phenomena that char- 
acterize star formation and feedback. So far, the coupling 
of a consistent modelling of the ISM with a reliable pre- 
diction of the physical properties of galactic winds turns 
out to be quite difficult. New generation numerical codes 
include often an explicit recipe for starting galactic winds 
in suitable situations. In these cases of defaillance, simu- 
lations have to rely on phenomenological prescriptions de- 
duced from theoretical models or inferred from the analy- 
sis of the local Universe. To date, the number of observed 
galaxies with on-going outflows is quite consistent with the 
theoretical predictions, but the correlation between the out- 
flow and galaxy properties i s difficul t to b e interpreted and 
shows no substantial trends (|Martinl dl999h jHeckman et alJ 
(2000 )). Recently, observations made by Martin 1 20051 ) and 
iMurrav. Quatert fc Thompson! l|2005l ) showed a good agree- 
ment with theoretical models for momentum-driven winds. 
In this model, radiation impinges on dust grains, transfer- 
ring this way its momentum. The latter is moving out of 
the galaxies and then may couple with the surrounding ga s 
and start an outflow dMurrav. Quatert fc Thompsonll2005r ). 
Based on this model, lOppenheimer fc Pavel (|2006h have run 
a set of numer ical simulations using the code GADGET2 
jSpringell [2005) with a span of feedback prescriptions for 
galactic winds. They demonstrated that momentum-driven 
outflows, calibrated via the local universe observations, are 
able to reproduce a wide variety of Civ absorbers statistics 
as observed in QSOs spectra. 

In this framework of these recent studies we have per- 
formed our study. We have used an output from a GAD- 
GET2 cosmological simulation that uses state-of-art treat- 
ment of multiphase ISM and winds, and keeps track of the 
contribution of Pop III stars to the overall metal budget. 
These simulations have been then analyzed with the inten- 
tion of providing quantitative and clear-cut diagnostics for 
observational investigation of outflow phenomena around 
high redshift galaxies; to this aim, and also to clarify the 
association between absorbing lines and galactic environ- 
ments, we have performed a spectroscopic analysis around a 
prototypical galaxy characterized by ongoing wind-activity. 

This paper is organized as follows: in §2, we briefly de- 
scribe the numerical simulation; we then move to the galaxy 
selection criteria and to the analysis of its gaseous environ- 
ments. In §3 we analyze the metal absorption features of the 
"wind region" observed in our spectra and we discuss briefly 
of the Lya forest around our target galaxy. In Section §4 we 
give then our conclusions. In addition, we provide in Ap- 
pendix A details about the method used to generate our 
spectra and a summary about the considered metal transi- 
tions. Finally, in Appendix B we show our synthesized spec- 
tra. 



2 NUMERICAL SIMULATIONS 

Our analysis is based on a hydrodynamic cosmological simu- 
lation. The simulation adopts a ACDM cosmological model 
with n M =0.3, fl A =0.7, ^(,=0.04, h = 0.7, and erg = 0.8; 
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periodic boundary conditions are adopted. The numerical 
code used for the simulation is a modifi ed version of GAD- 
G ET dSprineel. Yoshida fe WhitdkoOll l modified according 
to iMarri fc White! (|2003f T where the interested reader can 
find the details of the new algorithms introduced. In brief, 
additions concern: (i) an improved treatment of the mul- 
tiphase nature of the gas, and (ii) a more refined recipe 
for stellar feedback. The ISM is modeled as a two-phase 
medium, made of a cold, dense component and a warm, dif- 
fuse one. Once star formation occurs, the feedback energy is 
distributed to the surrounding particles with a fraction e c 
(e w ) incorporated in the cold (warm) phase. Values for e c 
and e w could be fixed from a complete theory of interstellar 
medium, or by specific numerical simulations; in our case 
these constants were both fixed to 0.4. The remaining 20% 
of the energy is assumed to be radiated away. 

The scheme also tracks the gas metal enrichment due to 
supernova (SN) explosions by considering a heavy element 
yield yz such that: 

5M Z = yzSM* (1) 

where SMz is the total metal mass produced per stellar mass 
(SM») formed. Once metals are produced, a fixed fraction 
/z=0.3 is distributed among the same hot neighbors which 
receive the feedback energy. Tracking the metallicity allows 
us to apply the so-called 'critical metallicity' criterion intro- 
duced by Schneider et al. (2002). According to this criterion, 
the Initial Mass Function (IMF) of star forming sites is bi- 
ased towards very massive (PopIII) stars if the metallicity 
is below Z cr = 10~ 4 Zq; above this threshold a switch to 
a standard stellar population (PopII) with a Salpeter IMF 
occurs. These stellar populations have distinct feedback pre- 
scriptions, depending on the explosion energies, e sn , and 
metal production efficiencies, yz of PopIII and PopII stars. 
Therefore, we adopt the following set of parameters: 

e*J? = 2.3 x 10 50 erg M _1 yi 11 = 0.637 (2) 

for PopIII stars, and 

e{' n = 4.0 x 10 48 erg Mq" 1 y' z ' = 0.004 (3) 

for PopII stars. These above parameters are not yet fully 
fixed by theory, as the outcome of stellar explosions depend 
e.g. on progenitor composition, nuclear reactions rates, ex- 
plosion energy. 

We have run simulations at two different resolutions. 
The low resolution simulation, used mostly as a control 
run, consists of 64 3 particles simulating a 7.0 h~ x comoving 
Mpc; the high resolution run is made up of 128 3 particles 
in a 10.5 ft -1 comoving Mpc. Both have a mass resolution 
of M ~ 10 6 Afo/i -1 . In both runs the evolution was fol- 
lowed including the effect of a UV-background (produced 
by QSOs and galaxies and subsequently filtered through 
the IGM) whose shape and amplitude were taken from 
lHaardt fc Madaul £l996). Unless otherwise specified we re- 
port results based on the high resolution run. 

2.1 Target galaxy 

As general test for our metal feedback recipe, we have ana- 
lyzed the low resolution run for the average metallicity of the 
IGM as traced by carbon, oxygen and silicon. In Fig. [T]wc 
show the relation of these abundances, relative to solar, as a 
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-1.99 


-2.00 4 



1 Schave et al. (2003) . 2 Simcoe et al. (20041 
3 Telfer et al. (20021 . 4 Aguirre et al. (2004) 

Table 1. Abundances values relative to solar from our low res- 
olution run: shown are the considered element (first column), 
transition used as tracer (second), median value retrieved from 
our simulation (third), observed values (forth). For a meaningful 
comparison with observations, we have included in our analysis 
only pixels with Lya and metal optical depth higher than 10~ 3 . 
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Figure 1. Abundances relative to solar for C (dotted), 
(dashed) and Si (dot-dashed) as a function of overdensity. The 
thick red curve represents the values for the total metallicity rel- 
ati ve to solar while the straight blue line show the relation found 
bv lSchave et al ] j2003h for [C/H] (eqHJl. 



function of the overdensity (8 = pj < p >). In order to com- 
pare with observations, we have limited our analysis to those 
pixels characterized by optical depths for both Lya and the 
metal transition > 10~ 3 . From the plot we deduce that from 
our simulation a small raise in the metal abundances should 
be expected at low overdensities where wind effects are more 
evident. The combined action of our thermal and chemical 
feedback seems to suggest also the fact that the [C/H] ratio 
could be the closest numerically tracer of the total metal 
content of the gas. In t he Fi gure we also show the relation 
found bv lSchave et al.l l|2003h for [C/H] vs. overdensity: 

[C/H] = -3.47±°'°6 + 0.65i° ;i4 x (log 5 - 0.5). (4) 

Tab.[T]compares the simulated mean abundances of the three 
elements/tracers with the observational values. From the 
joint analysis of Fig. [T] and Tab. [1] we conclude that our 
model is quite successful in reproducing both Schaye rela- 
tion and the level of pollution in the IGM as derived by 
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QSOs absorption lines. This strengthens our confidence on 
the implemented feedback scheme. 

In order to quantify and diagnose the effects of feedback 
and winds on galactic environments, the first step is to iden- 
tify a suitable representative target galaxy in the simulation 
box showing a strong wind. To this aim, we have searched in 
the high-resolution run output at z = 3.26 for high density 
(n > 10~ 2 cm -3 ), high temperature (T > 10 6 K) parti- 
cles. The latter condition ensures that we are dealing with 
a galaxy that is surrounded by shocked gas; however, such 
high temperatures can also be reached during virialization 
and/or infall, so they do not necessarily imply that the ob- 
ject is experiencing a wind phase. We thus complement the 
high-temperature condition with that of a high peculiar par- 
ticle velocity, as infalling material is likely slowed down by 
the interaction with the halo gas. The total mass of the se- 
lected target galaxy which fulfills the above listed conditions 
is « lO n M0, i.e. a mass typical for Lyman Break Galax- 
ies (LBG). In the chosen volume, the average metallicity is 
« W~ 2 Z e . 

2.2 The target galaxy environment 

In Fig. [2] we show 2D maps of different quantities (tem- 
perature and velocity field, ram-pressure, total metallicity, 
Pop III stars contribution to the latter) in a plane across 
the target galaxy center. The signature of the wind is clearly 
seen as a region of enhanced temperature in the correspond- 
ing map, with a typical (physical) size ~ 0.4 Mpc. The shock 
caused by the outflow has heated up the gas to a tempera- 
ture ~ 10 6 K, but the outflowing gas is being slowed down by 
the infalling material. We refer to this region as the (wind) 
bubble. 

The interplay between the infalling and outflowing ma- 
terial shapes the bubble geometry, producing significant 
deviations from the commonly assumed spherical symme- 
try. Instead, we find that the shape of the bubble is bet- 
ter described by a pancake-like geometry with a radius-to- 
thickness ratio of ~ 2 : 1. The bubble extends almost or- 
thogonally to the filament that hosts the galaxy, with 4 or 
more lobes caused by the infalling material, which we de- 
note as (infalling) streams. An interface between the wind 
bubble and a stream is clearly identified by the low pressure 
arc (right top panel) to the left of the bubble; there, the 
outflowing material is actually stopped by the ram-pressure 
of the infalling gas causing the formation of a thin, dense 
shell. The analogous shell formed at the top lobe is barely 
recognizable as it has already started to recollapse onto the 
galaxy. To the right hand side of the bubble instead the col- 
lapse evolution is more advanced and the infalling material 
is almost freely penetrating into the hot bubble. As a result, 
the velocity field is very complex with regions of outflowing 
and infalling motions with peculiar velocities up to 250 km 
s _1 . Inside the galactic environment, the infalling material 
is re-directed by the galactic activity and trapped within 
the hot bubble. From the maps shown here we cannot iden- 
tify any outflow channel through which the gas could escape 
the bubble. Inside the bubble, the material now is heated to 
temperatures T > 10 6 K. 

The spatial distribution of the metallicity (lower panels) 
is very similar to the temperature, clearly indicating that the 
bubble has been created by supernova explosions within the 



target galaxy. However, the metallicity pattern correspond- 
ing purely to PopIII stars requires a previous star formation 
activity which pollutes the streams. This is deduced from 
PopIII-enriched cold clumps falling toward the galaxy seen 
in the upper portion of the map. Such clumps have been 
produced at earlier cosmic times by progenitors of the tar- 
get galaxy and later became part of the accretion process 
feeding it. They are pushed and compressed by the pristine 
material lying behind them and streaming from the filament 
into the bubble (see temperature map). These kind of con- 
densed clumps in the bubble had enough time to cool down 
to temperature T < Tw = 10 5 ' 5 K; throughout the follow- 
ing we identify these structures as cold clumps. Note that 
also the filament contains metals, the chemical composition 
being dominated by the first generation of stars. We find 
that PopIII contribution to the total metal content of the 
filament is > 30%. Metal absorption arising in filamentary 
gas therefore is important to study the nucleosynthesis of 
the first objects. 

In summary, the galaxy environment is characterized 
by two different environments: (i) the filament in which the 
galaxy is embedded, and (ii) the wind environment, in turn 
constituted by the wind bubble, infalling streams and cold 
clumps accreting onto the galaxy. 



3 IONIZATION ANALYSIS 

Having selected the target galaxy and characterized its en- 
vironment, we now investigate their properties by means 
of synthetic absorption spectra that can be directly com- 
pared with observed QSO absorption-line data. To this aim, 
we have analyzed ~ 1600 lines of sight (LOS) through the 
selected volum^E The ionization equilibrium, under the in- 
fluence of the UV background adopted in the simulation, 
has been calculated in each of the (~ 100) cells into which 
each LOS has been segmented, using the code CLOUDY9£0 
llFerland et al . 1998) . For the following analysis we have used 
ei ther the Pixel Opt i cal D epth (POD) technique de scribed 
in lCowie fc Songailal i|l99Sl ) and lAguirre et al.1 i|2002l) . or the 
automated Voigt- profile fitting using the code AUTOVP 
l|Dave et al.lll997T ). depending on the specific question ad- 
dressed. For further details on the simulated spectra we refer 
to Appendix A. 

As a first step to search for possible tracers of circum- 
galactic high-temperature regions, we have analyzed the gas 
temperature distribution shown in Fig. [3] (upper panel). 
The horizontal width of these lines denotes the tempera- 
ture range in which the transmitted flux, F = e - ' 1 "', due to 
a given ion opacity is F < 0.95, a level denoted also by the 
dashed line in the lower panel. The vertical position of the 
lines indicates the fraction of cosmic baryons polluted with 
that ion (we have defined a pollution threshold by requiring 
a ionization fraction > 10~ 2,5 for each species) integrated 
over the above temperature range. The distribution shows 
a peak at T ~ 2 x 10 4 K, accounting for about 20% of the 
baryonic fraction in the volume; an almost equal amount of 
matter resides at T ^ 10 5 K. Civ traces almost the same 
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Figure 2. Maps of various quantities across the target galaxy center. Upper left: temperature and peculiar velocity field (arrows; the 
largest one corresponds to 250 km s -1 ); Upper right: ram-pressure; Bottom left: total metallicity; Bottom right: PopIII-to-total metallicity 
ratio. Sizes are in physical units. 



amount (about 0.6%) of baryons as Siiv, but it is sensitive 
to a slightly higher temperature regime. A similar trend is 
found for oxygen ions with Ovi spanning a relatively narrow 
range around the peak temperature (« 3 x 10 5 K) for col- 
lisional ionization of Ovi, whereas Ovu extends for almost 
1.5 dex. We find that 0.1-0.2% of all the baryons should be 
already contaminated with these species at z ~ 3.26. 

The mean transmitted flux for these ions varies as a 
function of temperature. These profiles are useful to un- 
derstand the physical origin of a given ion. The double 
peak feature seen for Civ, Siiv and Ovu, shows the distinct 
effect of the radiative and collisional ionization processes, 
with the high temperature peak being associated with col- 
lisionally ionization and the other one (typically less pro- 
nounced) resulting from photoionization. Therefore, for any 
given species, the relative strength of these two features gives 
a measure of the importance of photo- and collisional ion- 
ization. In particular, Civ has a higher contribution from 
collisional ionization than Siiv. An important point is that 
both Ovi and, even more evidently, Ovu are produced by 
collisional ionization, and therefore they represent particu- 



larly useful diagnostics to study the hot intergalactic envi- 
ronments like wind cavities. 

In the following we characterize in detail the absorption 
patterns associated with the two types of target galaxy en- 
vironments introduced above, i.e. the wind and the filament 
environment. 



3.1 Wind environment tracers 

As already pointed out, inside the volume enclosed by the 
supernova-driven shock three different gas phases/structures 
coexist: infalling streams, cold clumps and hot, outflowing 
gas. As seen from Fig. 3] in which we show the spatial dis- 
tribution of the ionization fractions of various ions, the bub- 
ble interior is characterized by a very low Hi abundance 
(log /hi < —6), caused by the enhanced collisional ioniza- 
tion rate associated with the high gas temperature. Such 
low abundance in the central regions is common to most 
of the species shown, with the only significant exception 
of Ovi. The latter species appears to be almost homoge- 
neously distributed in the volume, although its abundance 
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Figure 4. Ionization fractions (i.e, the fraction of the X element in the ionization state i with respect to its total abundance, fxi = 
nx i/nx) rnaps across the target galaxy center for various ions are shown. The dotted lines indicate the analyzed line of sight (LOS), 
numbered from 1 to 18 as labeled to the right of each panel. 
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Figure 3. Upper panel: temperature distribution of the baryonic 
fraction. Horizontal lines mark the temperature range in which 
the indicated ion produces a mean transmission of ~0.95; their 
vertical positions refer to the fraction of cosmic baryons polluted 
with that ion. Lower panel: mean transmitted flux, F = e~^ T \ 
for the different ionic species as a function of gas temperature. 



drops in isolated regions: the most noticeable example are 
the Ovi "voids" seen on the upper-left part of the distribu- 
tion map for this ion. Guided by the previous discussion on 
the nature of the different structures in the target galaxy en- 
vironment, we interpret such Ovi voids as cold clumps that 
condensed out of the expelled material and now falling back 
towards the galaxy. In addition to Ovi (which appears to be 
the best wind tracer) the bubble region is characterized by 
considerable abundances of Ovn and Ovm (for a detailed 
view see Appendix C, Fig. IB6l and Fig. IB7[) . Although cur- 
rently unavailable in high-redshift observations these ions 
would represent very reliable tracers of active outflow re- 
gions. 

In addition to hot gas (which has a very com- 
plex/turbulent velocity field, resulting in either outflowing 
or infalling flows) we have noticed the presence of cooler 
gas, either in form of infalling streams or isolated clumps. 
These cool phases are best traced by low-ionization species 
such as Cn, Cm, Sin, and Sim. For example, these ions have 
enhanced abundances in the cold gas clump previously iden- 
tified in the upper left region of the maps. Intermediate ion- 
ization ions (Civ and Siiv) are instead typically found in 
the interfaces between the cool and the hot gas phase. It is 
worth noticing that while Siiv has a more patchy distribu- 
tion concentrated around the clumps, Civ extends deeper 
into the bubble region and thus into higher temperature re- 
gions, as noticed previously (Fig. [S] left panel). A first qual- 
itative conclusion drawn from the maps is that the bubbles 
could be identified by searching for absorbers that simulta- 



neously have low Siiv/Civ (log(/suv//crv) < —2) and very 
high Ovi/Civ (log (/ovi //civ) > 1) ratios. These conclu- 
sions are visually supported by the two maps showing the 
spatial distribution of these ratios reported in Fig. [5] typical 
examples for this behavior are seen in the spectra obtained 
from LOS #10 and #14 shown in Fig. |S] where we present 
a selection of absorption spectra through the bubble envi- 
ronment along the LOS defined in Fig. [4] Both LOS exhibit 
strong Ovi and weak Civ absorption systems caused by in- 
falling (v > 0) material, but no Siiv is detected. Despite 
the higher Civ ionization fraction along LOS #13 (Fig. [4| 
no Civ absorption is seen along this sightline (Fig. [5} due 
to the very low gas densities sampled by that LOS. Stated 
differently, a non-detection of Civ absorption does not nec- 
essarily imply a deficiency of that ion. Instead, its distri- 
bution in physical and/or velocity space could be such that 
the corresponding absorption is smeared out. Only sightlines 
piercing regions of enhanced density (e.g. LOS #14) and/or 
located nearby the target galaxy (LOS #10) exhibit clear 
signs of Civ absorption. 

The cold clump itself instead is characterized by the 
presence of Sin, Sim and Siiv. It appears that the volume 
occupied by a given element (here Si) increases with its ion- 
ization state; the same conclusion holds for carbon ions Cn, 
Cm and Civ. Typical absorption patterns arising from the 
clump structures can be seen in LOS #15, #16 and #17 
shown in Fig. [U Absorption features arising from ions with 
similar ionization states (e.g., Sim and Cm) are well aligned 
in radial-velocity space, although some ions may be more 
spatially extended than others (e.g., Civ vs. Siiv). The con- 
tiguity with the hot gas in which the clump is immersed gives 
rise to a broad interface region and causes a Ovi absorption 
velocity shift of Av ~ 25 km s _1 with respect to Sin (LOS 
#16). It is worth noticing that along LOS #16, piercing the 
clump, doubly-ionized species (Cm and Sim) are heavily 
saturated (see Fig. |SJ). 

3.2 Comparison with the filament environment 

As we have discussed, the target galaxy and its wind bub- 
ble are embedded in a cosmic filament essentially aligned 
perpendicular to the plane of the maps shown in Fig. 1, i.e. 
its symmetry axis is parallel to the rc-direction. In order to 
investigate the spectral differences between the wind and 
the filament environments, we have compared LOS passing 
through the wind bubble and those intersecting the filament 
(hence at an x elevation larger than the bubble radius « 0.2 
physical Mpc). 

To anticipate the main point of the following analysis, 
we find very difficult to discriminate cosmic web structures 
from those created by the galactic wind purely from avail- 
able absorption data. The main reason is that cold clumps 
resulting from the fragmentation of the supernova-driven 
shell have properties that hardly differ from those of the 
filament hosting the system. 

To further elucidate this point, let us consider a LOS 
passing through the bubble on the y — z plane identified by 
x = 0; such line will cross the cool clumps/shell delimit- 
ing the hot cavity. Most of the absorption systems will arise 
from these dense structures rather than from the hot, rar- 
efied interior. Thus, although the latter gas has well defined 
abundance ratios (e.g. high Ovi/Civ, low Siiv/Civ as found 
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Figure 5. Spatial distribution of Siiv/Civ (left panel) and Ovi/Civ (right) ratios in the simulated volume around the target galaxy at 
z = 3.26. Sizes are in physical units. 



in the previous analysis of this phase) , their absorption sig- 
nature is too weak to be detected; on the contrary, the most 
evident absorption features come from the cool component. 
These will not be different from those obtained from the 
analysis of a LOS passing at higher x through the filament 
only, given that the physical conditions in the two structures 
are similar. 

It has to be noted that these similar physical conditions 
originate from different thermal histories. The cool clumps 
and shell are the product of the radiative cooling of the 
shocked bubble gas from high temperatures, T w 10 6 K. 
In the simulation accretion shocks are clearly s een in the 
filame nts, too (this fact was first pointed out bv lCen et al.l 
l| 19941 )) and these induce complex interface structures which 
harbor Ovi and other highly ionized species. Stated differ- 
ently, the effect of supernova-induced shock is not easily dis- 
entangled from shocks associated with the gravitational col- 
lapse of small scale structures giving rise to the overdensities 
we identify with the Lya fores10. 

A more quantitative conclusion on this issue might come 
from the POD analysis of the synthetic spectra presented in 
Fig. [7] Shown are the optical depths of various low- and 
high-ionization species plotted against the Hi Lya optical 
depth in the filament and in the wind environment. To cor- 
rect for saturation, we have complemented the Lya line with 
the weaker Ly/3 and hyy transitions. The two vertical dot- 
ted lines indicate the saturation point (i.e. when absorption 
reaches 85% of the continuum level) for Lya and Ly/3, re- 
spectively, showing from left to right the regions where these 
three transitions were used. To compare our results with ob- 
servations, we applied the same analysis to a sample of 28 
absorption systems found in 13 spectra from the VLT/UVES 
Large Program " The cosmic evolution of the IGM^ (J. Berg- 
eron et al.) which provides high S/N, high spectral resolution 
IGM data for z ~ 3 (see Tables ES [gp . 

Despite the fact that we are tracing high density re- 
gions, metals are detected down to Lya optical depth 



A similar effect has been noted by Kawata & Rauch (2007) 



~ —2.0. This is mostly due to the high signal-to- noise ra- 
tio [S/N(Crv)~ 150, S/N(Ovi)~ 98] at the redshift consid- 
ered. Unfortunately, a clear distinction in the behavior of 
the various ions in the two environments does not emerge: 
both the filament and the wind LOS opacities share al- 
most the same trend with the Lya optical depth. As ex- 
pected, Sin and Ovi are complementary tracers for the 
Hi absorption. Sin traces mostly the high-density clouds 
while Ovi is present in lower density regions. In simulated 
spectra half of the Hi absorbers, almost regardless of their 
physical properties, do not exhibit Ovi absorption. The ex- 
istence of these sub-sample lowers the mean Ovi optical 
depth and causes the second plateau at high Lya optical 
depth, where the number of high Ovi optical depths points 
has a drop. There is a good agreement between our simu- 
lated data and the observations for almost all ions, with the 
only exception of Siiv and, to a lesser extent, of Ovi. This 
is most probably due to a selection effect in our sample of 
absorption systems. For the Ovi our sample does not repre- 
sent properly the Ovi-free absorbers at high Lya opacities. 
As for Siiv, the shift of +0.75 dex seen between predicted 
opacities and the data is due to the fact that observed sys- 
tems are more biased towards intermediate column density 
absorbers and they are characterized by very shallow lines; 
while in our spectra almost all the Hi absorbers with column 
density log Nm — 13.8 are associated with relatively strong 
Siiv absorbers. 



In addition to the optical depth, the shape and width of 
the absorption systems provide additional physical informa- 
tion that might be in principle used to distinguish between 
absorbers in the bubble and in the filament. In fact, these 
environments are expected to be dynamically very different. 
In Fig. [5] we show the Doppler parameter distributions for 
the transitions considered (as obtained from Voigt profile fit- 
ting) in the filament and bubble environment. Once again, 
and for the reasons explained above, no significant differ- 
ences between the two distributions emerge for the consid- 
ered ions. 
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Figure 6. Absorption systems for the LOS #7-18 (labeled at the top of each panel) through the wind bubble; see Fig. [4] for their 
positions and ionization conditions. Negative (positive) velocities denote approaching (receding) flows with respect to the observer. 
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Figure 7. Pixel optical depth analysis results for the synthetic spectra. Plotted are the optical depths of various low- and high-ionization 
species against the Hi Lyo optical depth in the filament (solid lines) and in the wind (dashed-dot) environment; the starred line represents 
observational data from the VLT/UVES LP (see Tab, I All and Tab, IA2l , The two vertical lines mark the points at which Lya (leftmost) 
or Ly/3 (rightmost) absorption reaches 85% of the continuum level. Only bins with more than 50 data points are shown. 



3.2.1 Alternative strategies 

As pointed out before, the bubble environment is charac- 
terized by logT > 5.5 and highly ionized species such as 
Ovi, Ovu and Ovm (Fig. [2] [3}. Among these ions, only 
the Ovi transition falls within the UVES wavelength range, 
while those of Ovu (A = 21.6 A) and Ovm (18.97 A) are 
located in the X-ray band and will be observed by future 
X-ray facilities. It is interesting to briefly outline the gen- 
eral trends of the optical depth distribution also for these 
species. 

The optical depth for each transition has been 
computed using t he following in situ approximation 
|Theuns et alj|l998j ): 

r(fc) ~ n x (i)A px {a ° fXo) I <f G auss(k), (5) 

where the symbols are defined after Eq. lA2l and the physical 
properties of the gas are taken from the output of the numer- 
ical simulation. To check the sensitivity of the optical depth 
to temperature, as a test we have fixed the density to the 
typical values found in the bubble and in cool/filament re- 



gions, that is log(n) m —3.7 and log(n) as —2.7, respectively. 
In Fig. [9] we show the predicted Ovii/Ovin optical depths 
(normalized to the corresponding Ovi one) as a function of 
temperature for both environments. The optical depths for 
Ovu and Ovi become equal for temperatures logT « 5.6 in 
both environments; the ratio between Ovm and Ovi opti- 
cal depths becomes equal to unity at higher temperatures 
in the bubble region (log T ~ 6.0) than in the filament 
(logT ~ 5.7). Both Ovm and Ovu are then expected to 
provide larger opacities at high temperatures (logT ^ 5.5) 
with respect to Ovi, implying that these species would be 
particularly well suited to study the hot component of high 
redshift galactic winds. 

3.3 Additional filament properties 

Finally, we briefly present a statistical analysis of the Hi Lya 
forest in the filament region in which the target galaxy and 
its wind bubble are embedded. As the cosmic volume we are 
considering is biased by the presence of the relatively mas- 
sive target galaxy, we do not expect it to be representative of 
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Figure 8. Simulated Doppler parameter distribution for the different ions. Dotted (solid) lines refer to the filament (bubble) environment. 



the typical low-density environments sampled by Lya forest 
observations. Nevertheless, it is intriguing and instructive to 
evaluate the possible differences. 

To compare our sample of Hi absorbers with the ob- 
served properties of the Lya forest in high-redshi f t spec - 
tra we have used the data published bv iKim et all (|2002h . 
These consist of Lya absorbers analyzed by Voigt profile 
fitting of the spectra of 8 QSOs covering the Lya forest at 
redshift 1.5 < z < 3.6 obtained from VLT/UVES observa- 
tions. For consistency, we have only used a sub-sample of 
quasars with z > 3. The results of this comparison in terms 
of the neutral hydrogen column density, JVhi, and Doppler 
parameter distributions are shown in Figs. 1 101 111 The agree- 
ment in terms of the JVhi distribution is generally good in 
the range 13 < log JVhi < 15. The very low column density 
deficiency can be possibly caused by the different behav- 
ior of_th^_dat a anal ysis software used for the two studies. 
In IKim et all (|2002h the study was carried using VPFIT 
(Carswell et al.) compared to which AUTO VP seems to un- 



derestimate the number of low column density absorbers. 
Although data incompleteness can also play a role. As for 
the Doppler parameter, the mean of the simulated sample is 
shifted towards values that are lower than the observed one. 
This is plausible if we consider that in our volume we have 
large amounts of dense and cold gas (i.e. filament gas) caus- 
ing narrow absorption features due to the reduced thermal 
line broadening. 



4 CONCLUSIONS AND DISCUSSION 

In this work we have presented an analysis of artificial 
absorption- line spectra from the circumgalactic and inter- 
galactic gas of a wind-blowing galaxy selected from az~3 
output from a SPH simulation. The simulation includes both 
a multiphase treatment of the interstellar medium (ISM) 
and a con sistent picture for th e chemical enrichment of the 
Universe l|Bruscoli et al.ll2003l ). With our feedback scheme 
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Figure 9. Ratios between Ovu and Ovi (upper row), and Oviii and Ovi(lower row) optical depths. The left (right) column refers to 
the filament (bubble) region. The vertical red lines indicate the temperature where the two optical depth become equal. 
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Figure 10. Comparison between the Hi column density distri- 
bution of the absorbers in the simulated galaxy environment (red 
line) and the obs erved sample (diamonds). Data are taken from 
iKim et al] (|2002h . 



Figure 11. Comparison between the Hi Doppler parameter dis- 
tribution of the absorbers in the simulated galaxy environment 
(red line) and the ob served sample (di amonds). The observational 
data are taken from lKim et al.l l|2002l l. 



we are able to reproduce the observed values for [C/H], 
[0/H] and [Si/H] without requiring an overall constant pol- 
lution of the IGM as well as the relation bet ween [C/H] 
versus overdensity found bv lSchave et al.l (|2003l ) . 

The simulated object is comparable in mass with a LBG 
at this redshift. The spectra were generated to mimic the 
behavior of the UVES spectrograph. 

For the wind bubble we find that the interplay be- 
tween the cooler infalling material and the outflowing hot 



gas shapes a bubble-wind geometry, causing a topology 
more similar to a pancake like geometry developing orthog- 
onally to the filament hosting the galaxy, leaving the fila- 
ment almost undisturbed. We identify more than four lobes 
caused by the cold-bubble material, separated by the so- 
called infalling-channels. If this is typical for wind-blowing 
galaxies at high redshift, then the gaseous environment is 
determined by: (i) the cold filament hosting the galaxy, (ii) 
a mixture of hot gas (T > 10 5 K) filling up the cavity created 
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by the galactic wind, (iii) the almost pristine IGM infalling 
towards the galaxy and opposing to the hot-bubble expan- 
sion, (iv) the colder blown out material that could have had 
the time to fragment into smaller clumps. 

^From the spectral analysis we find that cold- 
bubble clumps are characterized by Sin, Sim, Cn and 
Cm absorption. High oxygen ions lines are not expected 
in these clumps, but the cooler gas (as indicated by 
Sin absorption) is enveloped by an Ovi phase. 

The outflowing gas from the galaxy is characterized by 
the high ionization states of oxygen, i.e., Ovi, Ovu, and 
Ovm. While only Ovi is observable in UV/optical absorp- 
tion, we show that the Ovu and the Ovm phase should give 
rise to absorption with substantial optical depths in the X- 
ray band. 

A physical environment similar to that of the cold 
clumps within the bubble is present in cooler regions in the 
surrounding intergalactic filament (i). Here, the cold and 
warm filamentary material produces low- and high ion ab- 
sorption patterns that mimic those found in the bubble envi- 
ronment. From an analysis of the optical depth and velocity 
width distribution of the absorption in both environments 
we find that there is no clear way to distinguish between 
these two spatially very distinct environments based on their 
absorption- line patterns. 

The comparison of our optical depths with the observed 
data shows a good agreement for almost all ions except for 
Ovi and Siiv. This is most likely due to the composition 
of our observed sample. In our simulated spectra, half of 
the LOSs showing strong Lya optical depths do not show 
appreciable Ovi absorption; while our observed sample lacks 
in these systems causing a higher mean Ovi optical depth. 
As concerning Siiv, our observed sample is mostly composed 
by intermediate column density absorbers. In our simulated 
spectra instead, higher Hi column densities absorbers are 
preferably associated with medium/high column densities 
Siiv absorbers. 

Our study suggests that it is quite difficult to discrim- 
inate between the various circumgalactic and intergalactic 
gaseous phases of a starbursting galaxy environment based 
solely on the absorption line characteristics. Enlarging the 
sample of target objects would probably increase the statis- 
tical weight of our conclusions, but we do expect that the 
general picture is preserved. The incidence of cold clumps 
on the spectra is expected to increase in more complex en- 
vironment, as for example if more filaments are intersecting 
each other. Along with a different metal distribution, bub- 
bles will be instead characterized by almost the same range 
of temperatures and densities; they will then give raise to 
very similar absorption patterns. Metal cooling can then be 
the most important missing ingredient of the present study. 
Non-equilibrium ionization effects could also leave an im- 
print on the cooled gas, which might still keep memory of 
the past hot and rarefied condition. We are planning to ex- 
plore these variables on future works. 

Fina lly, a brief comparison w ith tw o recently appeared 
studie s (|Oppenheimer fc Pavel l|2006l ); iKawata fc Rauchl 
l|2007l ) ) must be given. 

In lOppenheimer fc Pavel (|2006T l the authors mainly use 
statistics related to Civ absorption to compare their differ- 
ent feedback recipes with observational result. Consistently 



with their findings, we confirm that most of the Civ is of 
collisional origin; however, our results also r e produ ce the 
[C/H] vs. overdensity relation bv lSchave et all (120031 ). Most 
probably, th is discrepancy is due to a di fferent metallicity 
distribution. lOppenheimer fc Pavel l|2006h predict a metal- 
licity ~ 1 dex higher than our values; also their fiducial 
models show an almost flat pollution at high overdensities 
(see their Fig. 10) contrary to our continuously rising level. 
Still, both works predict that winds expand preferably in 
low density regions, leaving the filaments hosting the galaxy 

unaffected. 

As for IKawata fc Rauchl (|2007f ). their main finding is 
that Ovi is the ion most sensible feedback tracer. Their 
assumption of a constant - and somewhat unphysical - 
IGM pollution does not allow a direct comparison with our 
metallicity distribution. Nevertheless, by looking at their 
maps it seems that the bulk of the Ovi absorption arises 
from shocked filament material as well as from clumps con- 
densed from blown out gas, a scenario in agreement with 
the one proposed here. Repeating their search for strong 
Ovi absorbers (with flux > 0.8) associated with unsaturated 
Hi ones, we find that these systems are almost evenly pro- 
duced in the bubble environment and the filament one. We 
cannot then confirm the exclusive association of these sys- 
tems to winds. 
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APPENDIX A: SYNTHETIC SPECTRA 
Al Smoothing 

To calculate a spectrum along a given line of sight (LOS) 
through the simulation box it is necessary to recover the 
physical properties of the gas at each point along the path, 



0.10 




100 200 300 400 500 600 
n. part. 

Figure Al. Smoothing convergency test using a gaussian kernel 
with a FWHM a = 8.0 h" 1 kpc. 



by smoothing the various fields. For this, the line of sight has 
been divided into cells, to which the values of the gas den- 
sity, temperature, velocity, and metallicity are assigned by 
a weighted mean on the surrounding particles. The weight 
of each particle is determined by its distance from the LOS- 
cell using a specified kernel function. For consistency with 
the adopted SPH simulation, we have chosen a gaussian ker- 
nel with FWHM a = 8.0 h _1 kpc comoving. The number 
of particles used to derive the mean is fixed by the balance 
between convergence and computational expense. To deter- 
mine such number, we have recovered the various physical 
quantities for 600 cells along the LOS using a variable num- 
ber of particles (up to 600, the reference case). In Fig. IA1I 
we plot the results of such test in terms of the quantity 



t _ I \Xj ~ X 600 I \ 

w — \ j li 

"^600 



(Al) 



where Xf indicates the j'-th physical quantity recovered us- 
ing i particles, and the mean is calculated over the cells 
considered. As we see, an accuracy of < 6% is obtained by 
using 100 particles; by further taking into account the com- 
putational expense, we have fixed this value as the optimal 
one and have adopted it in the following. 

To further test the smoothing algorithm, we have recov- 
ered the physical properties of 1000 gas particles. Fig. IA2I 
shows a comparison between recovered and simulated val- 
ues. The accuracy of recovered quantities is reasonably good. 
The scatter in the distribution may be related to particles 
belonging to regions near the interfaces of virializing halos, 
the sharp transition across the shock front characterizing 
these regions been smoothed away by the recovery process. 
This interpretation is supported by the "wall" in the temper- 
ature distribution at T ~ W 4 K, i.e. the lowest temperature 
allowed by the adopted atomic cooling function. Moreover, 
the small filling factor of virialized regions tends to trans- 
late into a density overestimate. The velocity distribution 
seems instead more symmetric, as a result of the different 
kinetic behavior characterizing the collapsed haloes distri- 
bution. The metallicity field is particularly affected by the 
recovery process, as expected from the fact that no diffusion 
equations are solved for the heavy elements. As a result, pol- 
luted and pristine regions are separated by unphysical sharp 
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Figure A2. Smoothing recovery test for the density (left panel), temperature (middle) and velocity (right) fields. The dots represent 
the relation between the recovered and the simulated physical properties of 1000 particles in the high resolution simulation at redshift 
z ~ 3.0. The line represents the case of an ideal recovery algorithm. 



z 


QSO 




2.170 


Q0122 - 


380 


2.208 


PKS1448 


- 232 


2.280 


HE0001 - 


2340 


2.406 


Q0109 - 


3518 


2.406 


HE2217 - 


2818 


2.434 


HE1347 - 


2457 


2.740 


HE0151 - 


4326 


2.758 


Q0002 - 


422 



Table Al. QSOs spectra used to compute the models for S/N 
and S/E ratios for the UVES spectrograph. 



QSO 



2.134 
2.685 
2.661 
2.885 
3.054 



HE1341 - 1020 
PKS0329 - 255 

Q0453 - 423 
HE2347 - 4342 
HE0940 - 1050 



Table A2. QSOs spectra used to test the models for S/N and 
S/E ratios for the UVES spectrograph. 



boundaries. Results concerning the metallicity derived from 
our synthetic spectra should be taken with caution. 

A2 UVES models 

For our synthetic spectra we have chosen to mimic the be- 
havior of the UVES spectrograph mounted at the ESO Very 
Large Telescope (VLT). To obtain a reliable model for the 
signal to noise {S/N) and for the signal to error (S/E) ratio 
we have used a set of observations made as part of the Large 
Program " The cosmic evolution of the IGM" (J. Bergeron et 
al.) which provides high S/N, high spectral resolution IGM 
data for z ~ 3 (Tab. lATj) . 

Considering that the maximum restframe absorption 
wavelength used in our spectra is A ~ 1550 A, we have fixed 
the wavelength range for our models to 3300 < A < 7400 A. 



Hi 



I 

3000 



Figure A3. Fiducial model of the S/N for the UVES spectro- 
graph. In black we show the observed data, along with their vari- 
ance; while the red line refers to our polynomial fit. We can note 
the different density of data with A, this reflects the varying num- 
ber density of the absorbing clouds. 



A2.1 S/N model 

To obtain the S/N dependence with wavelength we have 
searched (by eye) for unabsorbed regions in the UVES spec- 
tra listed in Tab. IA1I The data obtained then were binned 
in 0.25 A bins, so that every bin is formed by 10 pixels at 
the UVES resolution of ~ 0.025 A per pixel. We then have 
considered the mean value as estimate for the signal from 
each pixel/bin, the noise was estimated by a simple subtrac- 
tion. The triplet of values (X-S-N) obtained in this way were 
collected in a single large sample sorted by wavelengths. To 
reduce the data used for the fitting, we have divided the fi- 
nal data set in bins of 100 values. For each bin we then have 
estimated the mean value and the scatter for the S/N by fit- 
ting the distribution of values using a gaussian. The values 
obtained were then fitted using a polynomial function (Fig. 
IA3[) . In Fig. IA4I we show the test of our fiducial model for 
the S/N using the QSOs given in Tab.|A2] We have plotted 
the distribution of the absolute scatter between predicted 
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Figure A5. Dependence of the S/E with wavelength. The black 
dots refer to the observed data, while the blue line represents our 
polynomial fit. Red dots represent the predicted ratios using the 
Poissonian approximation. 

and observed 5*/Ajf|. As we can see, also if the distribution 
reaches high values; the rapid decrease of the observed his- 
togram satisfies our model. 

A2.2 S/E model 

To obtain a model for the signal-to-error ratio, we have 
treated the observed data as described in the previous sec- 
tion. In Fig. IA5I we show the observed data (black dots) 
along with our model (blue dots) and a Poissonian model 
(red dots). We see how the latter approximation (usually 
considered) overestimates the errors associated with fluxes. 
^,From Fig. IA6I we can see that both the approximations 
overestimate the errors associated with fluxes. However, our 
model can be considered as a better low-limit approximation 
to the real ratio. 

4 The data were treated in the same way as for the model. 



A 2. 3 Binning issues 

To assess the dependencies of our S/N model, we have car- 
ried out the binning procedure using different bin sizes. In 
Fig. IA7I we show on the left our results for binning sizes of 
0.25, 0.50, 0.75, 1.0, 2.0, 4.0 A, respectively, as dots with 
colors from black to red. The right plot shows the mean rel- 
ative scatter of each data set, from the reference one (that 
relative to a 0.25 A binning). 

These two plots demonstrate that the main consequence 
of rebinning is a lowering of the observed S/N ratio, most 
probably due to the intrinsic continuum variation inside the 
bin. Since we do not have fully analyzed spectra, and con- 
sidering that every new set of S/N data lies into the error 
bars of our reference-data set, we can consider our model as 
acceptable. 

A3 Spectrum generation and degradation 

To obtain the absorption arising along the LOS, that means 
the optical depth r for each cell, we have used the Eq. IA2I 
|Theuns et alj|l998h : 

JVpx 

r(fc) = yV(i)Ap, ( 7f A ° } c <S>voi B t(k,i) (A2) 

i — 1 

n x {€) = (Z//(*)y//, e (<)+Zm(t)ym,«(*))n(i)/(A ie) x+)(A3) 
where: 

(i) k is the cell where we want to calculate the optical 
depth; 

(ii) the summation is over all the cells, N px , in each LOS; 

(iii) n x (i) is the numeric density in the cell i of the ele- 
ment x with ionization degree \+\ 

(iv) n(i) is the total gas density in the cell i; 

(v) Yn, x and Ym, x are the yields for the element x, re- 
spectively for PopII and PopIII stars (Tab. DO)) ; 

(vi) f(A x ,x+) is the ionization fraction; 

(vii) A px is the linear dimension of the cell; 

(viii) z is the redshift; 

(ix) ao, /, Ao are atomic constants, respectively: (To is 
linked to the Thomson cross section, <jt, as follow ao = 
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3000 4500 6000 7500 12 3 4 

A (A) AA bin (A) 

Figure A7. Dependencies of the S/N model on the binning size. On the left we show the observed S/N ratio for different wavelengths. 
Colors, from black to red, refer to a growing binning size (0.25, 0.50, 0.75, 1.0, 2.0, 4.0 A). For the first binning size (our working size) 
we show also the relative errors. The right plot shows the observed mean relative scatter of each set from the reference one as a function 
of the binning size. 



Element 


PopIII 


PopII 


C 


0.0488 


0.109 





0.504 


0.693 


Si 


0.211 


0.0626 



Table A3. Yields for t h e element included in our spectra 
llWooslev fc Weaverl Jl995h . iHeger fc Wooslevl ll2002h '). 



Vjf&T, and / is the oscillator strength for the transition 
(x) b x is the Doppler parameter for the element x defined 

as: 



b x (i) 



2kT(i) 



12.845 



km s 



(A4) 



where: k is the Boltzmann constant, T(i) is the temperature 
of the i-th cell in Kelvir0, m x and A x are respectively the 
mass and the atomic number of the element x; 

(xi) &voi g t(k,i) is the statistical weight that determine 
the contribution of the cell i-th to the absorption of the cell 
fc-th, in our case we have chosen to use a Voigt profile. 

The absorption wavelength associated to each r(fe) is 
obtained considering the rest-frame absorption wavelength 
at the redshift of the numerical outpufl and the total veloc- 
ity of the cell. To compare with the data, it is necessary to 
degrade the raw spectrum to account for finite instrumen- 
tal resolution and noise. To this aim we convolve the spec- 
tra with the (UVES) instrumental profile assumed to be a 
gaussian of width ~ 6.6 km s" 1 . The noise contribution to 
the flux in each bin is obtained by considering a gaussian 
noise with zero mean and with width determined by our 
S/N model. The same is done for the error associated with 
each flux value, fn Tables [A"4l and lA"5l we list the atomic pa- 
rameters for the atomic transition included in our simulated 



spectra. Here we show respectively in the first column the 
designation for the absorption followed by the rest frame ab- 
sorption wavelength, oscillator strength, gamma value and 
the second absorption wavelength in case of a doublet. 



APPENDIX B: SAMPLE SPECTRA 



In the following Figs. IBll - IBlOl we show the simulated spec- 
tra along a slice of the simulated volume containing the tar- 
get galaxy. 



5 The atomic para meters for each transition were taken from 
IVerner et aTl l ll994h . 

6 T 4 = T/10 4 K. 

7 The redshift difference between two cells is less than 0.01%. 
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T 

Ion 


Ai(A) 


f 
f 


„, /in9 
7/10 


A 2 (A) 


Civl548 


1548.195 


0.190 


0.2650 


1550.770 


Oil302 


1302.168 


0.0504 


0.5750 




Ovil031 


1031.926 


0.133 


0.4.125 


1037.617 


Siiil526 


1526.707 


0.132 


1.9600 




Siiil304 


1304.370 


0.0871 


1.7200 




Sinl260 


1260.422 


1.180 


2.5330 




Siiill93 


1193.290 


0.584 


3.4950 




Siiill90 


1190.416 


0.293 


3.5030 




Siml206 


1206.500 


0.0168 


2.550 




Siivl393 


1393.755 


0.524 


0.88250 


1402.770 



Table A4. Atomic parameters for the metals transitions included in our simulated spectra (tVerner et al. | [l994h . A2 refers to the second 
component in doublets. 



Ion 


Ai(A) 




f 






7/10 7 


HI1216 


1215.670 





416 






62.650 


HI1026 


1025.722 





0791 






18.970 


Hi972 


972.537 





0290 






8.126 


Hi950 


949.743 





0139 






7.640 


Hi938 


937.804 





780 ■ 


10" 


2 


4.4230 


Hi931 


930.748 





481 ■ 


10" 


2 


1.236 


Hi926 


926.226 





318 ■ 


10- 


2 


0.8249 


Hi923 


923.150 





222 ■ 


10- 


2 


0.5782 


Hi921 


920.963 





160 ■ 


10" 


2 


0.4208 


Hi919 


919.351 





120 ■ 


10" 


2 


0.3158 


Hi918 


918.129 





921 ■ 


10" 


3 


0.2430 


Hi917 


917.181 





723 ■ 


10- 


3 


0.1910 


Hi916 


916.429 





577 ■ 


10" 


3 


0.1529 



Table A5. Atomic parameters for the Hi included in our simulated spectra jVerner et al. Ill994h . 
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Figure Bl. The map shows for the ionization fraction of the ion indicated by the label. The dotted lines identify by progressive number 
the LOS whose spectra are shown in the upper panels, For ions not discussed in detail, we show in the upper panel the expected column 
density derived from an integration along each LOS. 
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Figure B2. 
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Figure B3. 
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Figure B4. 
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Figure B5. 
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Figure B6. 




Figure B7. 



A. P. M. Fangano et al. 




Figure B8 
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Figure B9. 
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